Robust chaos with variable Lyapunov exponent in smooth one-dimensional maps 
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We present several new easy ways of generating smooth one-dimensional maps displaying robust chaos, i.e., 
chaos for whole intervals of the parameter. Unlike what happens with previous methods, the Lyapunov expo- 
nent of the maps constructed here varies widely with the parameter. We show that the condition of negative 
Schwarzian derivative, which was used in previous works, is not a necessary condition for robust chaos. Fi- 
nally we show that the maps constructed in previous works have always the Lyapunov exponent In 2 because 
they are conjugated to each other and to the tent map by means of smooth homeomorphisms. In the methods 
presented here, the maps have variable Lyapunov coefficients because they are conjugated through non-smooth 
homeomorphisms similar to Minkowski's question mark function. 

PACS numbers: 05.45.Ac, 05.45.-a 



I. INTRODUCTION 



Many families of smooth maps display fragile chaos, which 
may be destroyed by arbitrarily small changes of the param- 
eter. For instance, the discrete dynamical system generated 
by the logistic map, Xn+i — I^Xn (1 — a;„), is chaotic for 
/i = 4, but the values of the parameter /i for which the at- 
tractor is periodic are dense in the interval [0,4] fl']. In con- 
sequence, when such a family is used to describe a phys- 
ical system, it may be impossible to decide on theoretical 
grounds whether the actual behavior of the system will be 
chaotic or periodic for some parameter value, which is nec- 
essarily known only approximately. Furthermore, some prac- 
tical applications, such as encoding messages [2], require re- 
liable chaotic behavior. 

Piecewise smooth maps may show robust chaos and they 
have been used to describe a circuit with robust chaotic out- 
put lUt]. Although for some time it was conjectured that one- 
dimensional maps should be piecewise smooth to display ro- 
bust chaos Ii3i,^4], Andrecut and Ali first found a smooth map 
lUt] and later a method of generating smooth maps |6] whose 
evolution is chaotic for whole intervals of the parameter. 

The purpose of this work is twofold: we want to ex- 
plore other easy ways of generating robust chaos in one- 
dimensional smooth maps and to check whether the condi- 
tion of negative Schwarzian derivative satisfied by the maps 
of Refs. Jslla] is a necessary one. Unlike in previous methods, 
the Lyapunov exponent of the maps explored here takes rather 
different values depending on the value of the parameter This 
property might be an advantage in some applications. We will 
see in Sect. |IV] that the reason of this dependence lies in the 
different ways in which maps generated by each method are 
conjugate to each other. 

We will consider maps on a finite interval [a, b], which for 
commodity will be reduced to [0, 1] by means of a linear trans- 
formation. 



II. ROBUST CHAOS WITH NEGATIVE SCHWARZIAN 
DERIVATIVE 

The Schwarzian derivative of the function / is defined as 



S/(x) 



(1) 



Since Singer used it in the study of the bifurcations of maps 
of the interval fT*], a key assumption in many theorems on the 
dynamics of one-dimensional discrete dynamical systems is 
that the Schwarzian derivative of the map is negative along 
the whole interval. 

In the following we will take advantage of the fact that 
the Schwarzian derivative is invariant under linear fractional 
transformations [8J, i.e., that for constants a, b, c and d one 
has 



afix) + b 
cf{x) + d 



= Sf{x) 



(2) 



Our starting point will be a map / : [0, 1] [0, 1] of class 
C^, which we assume to be 'S-unimodal,' i.e., which satisfies 
/(O) = /(I) = 0, has a single critical point at c e (0, 1), and 
negative Schwarzian derivative everywhere in [0,1]. Clearly / 
increases from its null value at a; = until it reaches its max- 
imum at a; = c, and then decreases until becoming again at 
X = 1. According to Singer's theorem [7, 8], the discrete dy- 
namical system a;„+i — f{xn) has at most one stable periodic 
orbit, which when exists attracts the critical point x = c. 

For any map / with the properties above mentioned we will 
construct the following one-parameter family of maps; 



frix) = 



{l + r)f{x) 
f{c)+rf{xy 



(-1 <r < oo). 



(3) 



(Notice that fr does not change if one multiplies / by any 
constant.) 

By using Q and 
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fr(.x) 



(l + r)/(c)f(x) 
[/(c) + r/(x)]2 



(4) 
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FIG. 1: Lyapunov exponent of the map ((Sjl for j[x) = x[\ — x) 



one can readily check that /r is also S-unimodal for all r > 
-1. 

Now, for 



r > To 



fjc) 
/'(O) 



1, 



(5) 



the origin is a unstable fixed point, because then (0) > 1 
and, furthermore, the dynamical system Xn+i = fr{xn) has 
no stable periodic orbit, because the critical point goes, in two 
steps, to the unstable origin: ff.[c) = 0. 

To check that the dynamical system is chaotic for all r > tq 
on can compute numerically the Lyapunov exponent 



TV 



^= J™ 5]ln|/;(x„)|. 

N—>-oc — ^ 



(6) 



With the maps generated from the logistic map /(x) — 
x(l — x) by means of (|3]l, one gets the values displayed in 
Fig.[T] They are negative for r < tq — —3/4, because then 
the origin is an attractor, and become positive at r > ro as 
the generic orbit wanders chaotically around the whole inter- 
val. In the maps generated in Refs. Jslla] the Lyapunov coeffi- 
cient was always In 2 or very close, according to the numerical 
simulations. In the maps presented above, the Lyapunov ex- 
ponent varies with r in a continuous way. If /(x) = a;(l — x) 
the maximum value of the Lyapunov exponent of /r(a;) is 
Amax = In 2. This value is reached at r = 0, which corre- 
sponds to the well known case Xn^\ — 4xn {1 — Xn), which 
in turn is conjugate to the tent map defined as f{x) — 2x for 
< X < 1/2 and f{x) = 2 - 2x for 1/2 < a; < 1 H. 
Similar graphs, with < A < Amax = In 2, are obtained, 
for instance, for the asymmetric map f{x) — x{l — x^) and 
for f{x) — sinTTx, although the bifurcation value tq and the 
location of the maximum are different. 

One can also compute numerically the natural invariant 
measure dfi = Pr{x) dx by using the Frobenius-Perron equa- 
tion iQl] satisfied by the natural invariant density Pr(x): 



Prix) 



Pr (yi) , Pr (y2) 



I/; (2/2)1' 



(7) 



FIG. 2: Natural invariant density of the map ([3} for f{x) — x{l — x) 
and r = (dashed line) and r — 1 (continuous line). 



and 2/2 are the preimages of x, i.e., fr (yi) = 
X. In the case of f{x) = x{l — x) it is well 



where yi 

fr{y2) = 

known that for the logistic map /o the natural invariant den- 

sity is po(2^) = [7r'^a;(l — x)] . It is displayed in Fig. |2] 
along with the natural invariant density for /i (x). Only in the 
first case (for r = 0) is the natural invariant density symmetric 
around the critical point x = 1/2. Similar results are obtained 
with other choices of f(x). 

When constructing smooth maps by using the method of 
Andrecut and Ali |6] or the one provided by Eq. ([3]), robust 
chaos is guaranteed by Singer's theorem; but they are by no 
means the only way to get chaos for an interval of the parame- 
ter For instance, we have been exploring the family generated 
from a S-unimodal map f{x) by the expression 



/(c) 



(8) 



The Schwarzian derivative S fr has a rather involved expres- 
sion which makes difficult, if not impossible, a general analy- 
sis. However, selecting the logistic map f{x) — x{l — x), it 
is easy to see that the corresponding fr{x) is S-unimodal for 
—4 < r < 4 and displays robust chaos for —3 < r < 4. The 
plot of the corresponding Lyapunov exponent is very simi- 
lar to that of Fig. [T] (including the location and the value of 
its maximum), except for the fact that the bifurcation hap- 
pens at r = —3. Similar results have been obtained with 
f{x) = X (l — x^) and f{x) = sinvrx. 



III. ROBUST CHAOS WITH POSITIVE SCHWARZIAN 
DERIVATIVE 

All the maps discussed above, as well as those of Andrecut 
and Ali |5, 6] and the 'B-ExponentiaF map of Ref. |9], sat- 
isfy the condition of negative Schwarzian derivative. This is a 
very powerful condition, but also rather restrictive and can be 
destroyed by a smooth change of the x coordinate [10] or by 
small perturbations [8]. In consequence, it may be of practi- 
cal interest to find robust chaos in one-dimensional maps even 
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FIG. 3; Lyapunov exponent of the map l|3ll for Singer's function 
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FIG. 4: Lyapunov exponent of the map l llOt for the logistic map. 

when that condition is not satisfied. We will see in the follow- 
ing that the condition is not necessary to have robust chaos in 
one-dimensional smooth maps. 

Let us first consider Singer's function 

fix) = 7.86a; - 23.31x2 + 28.75x^ - 13.3x^ (9) 

and the map fr generated from it by means of (O. Since 
f{x) has a positive Schwarzian derivative in a subinterval of 
[0, 1], exactly the same happens with fr{x). However, a nu- 
merical computation of the Lyapunov exponent shows that 
fr{x) has robust chaos after the origin becomes unstable at 
r « —0.88156. In Fig.[3]we can see that the maximum value 
of the Lyapunov coefficient is in this case Amax ~ 0.62, i.e., 
somewhat smaller than the maximum value In 2 obtained in 
all previous examples. As happens in those examples, there 
is no attractor as the generic orbit wanders chaotically around 
the whole interval. 

We have also explored the one-parameter family of maps 

/.(x)^(^My, (^>o), (10) 

for some choices of f{x). 



In the case of the logistic map f{x) = x{l — x), the map 
fr{x) is S-unimodal only when r — 1. For r > 1 the function 
has a minimum at the origin. In consequence, x = is a 
stable fixed point that attracts the generic orbit, after a chaotic 
transient, which may be very long for values of r just above 1, 
for then the basin of attraction of a; = is tiny. For < ?■ < 1 
the Schwarzian derivative is positive near the origin and x = 
1. For instance, 

S/r(x)^i^ asx^O. (11) 

Furthermore, the map is not even C"^ in that case, because its 
first derivative goes to infinity at a; = 0, 1. However, we can 
see in Fig.|4]that robust chaos arises after the fixed point lo- 
cated in the interval (1/2, 1) becomes unstable at r w 0.1759. 
The maximum Lyapunov exponent is again hi 2 and is reached 
atr — 1, when we recover the logistic map /(x) = 4a;(l — x). 
Similar results have been obtained with f{x) = a; (l — a;^) 
and f{x) = sinTTx. 

The examples discussed in this section suggest that ro- 
bust chaos may not be an unusual property of smooth one- 
dimensional maps, even when the condition of negative 
Schwarzian derivative is not satisfied. 



IV. ROBUST CHAOS AND CONJUGATE MAPS 

All the maps generated here and in previous works have 
qualitatively similar dynamics: the solution wanders around 
the whole interval in a chaotic way. Furthermore, the graphs 
of all the maps g — fr look rather similar: they start from 
.g(0) = 0, increase monotonically until g{c) = 1 and the de- 
crease monotonically until g{l) — 0. This suggest all the 
maps are conjugate 1 1] to each other, i.e., given two of these 
maps, g and g, there exist a homeomorphism on [0, 1] such 
that g — (/)o go(j>~^. In other words, there exists a continuous 
change of variables x x — (pix), with continuous inverse, 
such that 

^[g{x)]^g[^{x)], Va;e[0,l]. (12) 

The dynamical systems Xn+i = g (a;„) and a;„+i g (a;„) 
have essentially equivalent dynamics (for instance, if the un- 
stable periodic orbits are dense for g the same will happen for 
g). If, additionally, (j) and (f>~^ are smooth both maps have the 
same Lyapunov coefficient [fit]. 

This property provides a simple method of constructing 
families of maps with robust chaos and constant Lyapunov 
exponent: take a chaotic map f{x) and a smooth homeomor- 
phism (j>{x) on [0, 1] depending continuously on a parameter 
r. Them fr = 4>o f<j)^^ will have the same Lyapunov expo- 
nent for all values of r. For instance, using f{x) = 4a;(l — a:) 
and (j>{x) = x^ for r > we get 

fr{x)^4'-x[l^x^^-y , (r>0), (13) 

whose Lyapunov exponent will be A = In 2 for all r > 0. 
However, this is not S-unimodal, except for r — 1, because 
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FIG. 5: Change of variables x = <j){^) when g — 4x{l — x) and g is 
i fTel l with r = 10. 

its Schwarzian derivative becomes positive near a; = for 
r > 1 and near x = 1 for < r < 1. (Moreover, 
> for < r < 1.) We see again that the condi- 
tion of negative Schwarzian derivative is not conserved by 
smooth changes of coordinates and is not necessary for ro- 
bust chaos. Notice that, by using the well known results cor- 
responding to the full logistic map lit], we can write explic- 
itly the solution of the dynamical system driven by the map 

(fTSl l as Xn — sin^'^ ^2" arcsina;^^'^^ and its natural invari- 

ant density as p{x) = [n^r^x^ (x^^/^ — l)] . The 'B- 
Exponential' map of Ref. 19] also is conjugate to the logistic 
map g{x) — 4a;(l — x). 

On the other hand, given two maps, g and g, with the quali- 
tative properties mentioned at the start of this section, one can 
use a method of successive approximations to construct the 
function 0, if it exists. One may proceed as follows: 

M^) = 2;, (14) 
cj)n+i{x) = g-' (5^)] , (n = 0, 1, 2, . . .). (15) 

The inverse function g^^ is two-valued in this kind of map, 
but the right preimage is given by the condition that if c and c 
are the critical points of g and g then > c when x > c. 
The method can be checked by computing cj)^^ in the same 
way to make sure that g and cj) o g o cj)^^ agree to the desired 
accuracy. 

We have found that the method converges quickly when 

tx) = 4:x{l — x) and g is one of the maps generated in Refs. 
la]. For instance, if we choose g{x) as given by the map 

1 _ „-x(l-x) 

of Ref. |6], with r — 10 and g{x) = 4a;(l — x), we get the 
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FIG. 6: Change of variables x — ^(a::) when g = 4a;(l — x) and g is 
with r = and f{x) is given by 

result of Figure |5] Similar results are obtained with other val- 
ues of r and for the map in Ref. fs*]. All these maps are thus 
conjugate to f{x) = Ax{l — x) and, in consequence l^lj, to 
the tent map. Since the function </> and its inverse are smooth, 
all these maps share the Lyapunov exponent A — In 2. 

We have checked numerically that also the maps generated 
from f{x) = a;(l — a;) by means of the different methods pre- 
sented in this work are conjugate to f{x) = 4a;(l — x) and, 
thus, to the tent map. But there is a crucial difference: al- 
though (f) and (f)^ ^ are continuous they are not smooth enough 
for the two maps to share the same Lyapunov exponent. For 
instance, in Fig.|6]we have chosen g — 4a:(l — x) and g as 
given by (|3]l, with r = 0, for Singer's map (|9]l. It is clear there 
that 4>'{x) vanishes at some points and, since the graph of cf)^^ 
is obtained by exchanging the axes of Fig. |6] the derivative of 
4)^^{x) is infinite at those points. 

This explains why the corresponding Lyapunov exponents 
are different. From this point of view, one can understand the 
methods of previous works as an easy way to construct one- 
parameter families of conjugate maps by means of smooth 
homeomorphisms that guarantee the conservation of the Lya- 
punov exponent, while the methods discussed here are easy 
ways of constructing conjugate maps with Lyapunov expo- 
nents varying in a continuous way, since they are conjugated 
by non-smooth homeomorphisms. 

The complex structure of the graph of (j){x) in Fig.|6]can be 
explored by zooming in on small parts of it. For instance, in 
Fig- Hone can seen the function in the interval [0.5, 0.51]. (A 
similar zoom of Fig.|5]reveals a smooth structure.) 

In fact, the graph in Fig.|6] as well as the remaining graphs 
we have obtained with the maps generated by the methods 
presented in this work, looks very similar to the graph of 
Minkowski's question mark function li 1 IJ shown in Fig. |8] 
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Lyapunov exponents of the tent map (A = In 2) and the Farey 
map (A = 0). It also explains why the natural invariant density 
of the latter map is not normalizable: p{x) (x x^^. 

The numerical evidence we have found and the dependence 
on the parameter of the Lyapunov coefficient strongly support 
the conjecture that the functions (/> conjugating pairs of maps 
generated by the methods described in this work are not dif- 
ferentiable at an infinite number of points, probably almost 
everywhere. 

On the other hand, the fact that Minkowski's question mark 
function can be recursively constructed by using the Farey se- 
quence and continuity 1 15] suggests an alternative method to 
construct for functions g and g. One starts from the critical 
point .To = c, since we know j/o = ^(2^0) ~ 4>{c) = c. Then 
for each pair (x„,?/„ = (t){xn)) already computed, one can 
calculate two new pairs 



+ 1, Vn+l = (piXn+l)) = {g±^ {Xn) , ff±^ {y,i)) , (17) 



FIG. 7: Zoom of a 0.01 x 0.01 square of Fig.|6] 




FIG. 8: Minkowski's question mark function. 



The resemblance is even more striking when both g and g are 
symmetric around the point x = 1/2. 

Minkowski's ?{x) function is continuous with continuous 
inverse, strictly increasing and its derivative is zero almost ev- 
erywhere and infinite or undefined otherwise fTJ, 13]. The 
question mark function is the homeomorphism cj) conjugat- 
ing the tent map and the Farey map il4ll defined as g{x) = 
x/{l - x) for < x < 1/2 and g{x) = (1 - x)/x for 
1/2 < X < 1. We have used this fact to check the accuracy 
of the method of successive approximations given by (fT4]i- 
(fTsT l. The fact that ?(a;) is not smooth explains the different 



where gZ^{x) is the value y satisfying g{y) = x and y < c, 
while y = g^^{x) is given by the conditions g{y) — x and 
y > c. Analogous definitions are used for g^^. We have 
checked that applying recursively (fTTI i one obtains again Figs. 
|6]and|2] The method also works for other pairs of maps con- 
structed by means of Q, (|8j or (fTOl i. 



V. FINAL COMMENTS 

In previous examples — including those of Refs. IHI^ but 
excluding ( fT3] l — the maximum is located at the same point 
for all values of r; but this is not a necessary condition. Let us 
consider the one-parameter family of maps 



frix) 



~ /(c) 



{0<r< 1), 



(18) 



which is obtained from family ( fTOb by means of the smooth 
homeomorphism (/)(a;) — x^l"^ . For instance, if /(x) = a;(l — 
x\ the maximum of (fTSl l is located at x = 0(1/2) = 2^^/^ 
and the Lyapunov exponent is that of Fig. |4] and its maximum 
value In 2 is reached again when r = 1 and we recover the full 
logistic map j\{x) = 4a;(l — x). 

We have also considered the following family of maps: 



_ /[(l + r)x 



/(c) 



(-1 < r < 1) 



(19) 



If /(x) — x(\ — x), the maximum is located at x = 
( Vl + T^ — 1 + r) /(2r) and the origin becomes unstable for 
r ~ —3/4. Again the maximum Lyapunov coefficient is In 2, 
but it remains very close to this value for a large parameter 
interval, as shown in Fig.|9] 

In all the examples considered above, as well as in those of 
Refs. |5, 6], the Lyapunov exponent is never higher than In 2; 
but it is easy to get other maximum values by changing the 
starting map f{x). Let consider only a simple example. The 
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FIG. 9: Lyapunov exponent of the map il9\ for the logistic map. 
piecewise linear map 

{3a;, if0<x<l/3; 
2 -3a;; if 1/3 < a; < 2/3; (20) 
3a; -2, if2/3<a;<l 

has |5'(a;)| ~ 3, except at a; = 1/3, 2/3. In consequence, 
its Lyapunov exponent is In 3 and its natural invariant density 
p{x) = 1. If we use the change of variables x ~ (p{x) = 
sin^(7ra;/2), the conjugate map f = (f>o g o (p^^ is 

f{x)^xiAx-3)\ (21) 



The Lyapunov exponent of / is In 3, because is smooth. 
Since (p is precisely the map conjugating the tent map and 
the full logistic map g{x) — 4a;(l — a;), the natural invariant 
density of the later is also that of /. 

Although this map is qualitatively different from those con- 
sidered above (for instance, it has two critical points), the Lya- 
punov coefficient of the coiTesponding family (|3]l looks much 
like that in Fig. [T] except for the fact that the origin becomes 
unstable at r = —8/9 and that the maximum value at r = 
is now Amax = In 3. If one computes the homeomorphism 
conjugating two maps of the family by a trivial extension of 
the method (fT4b-(fT5]), one obtains a graph similar to that of 
Fig. |6l The same happens if one uses family ( fTOl ). in which 
case the Lyapunov coefficient is similar to that of Fig. 21 with 
the bifurcation at r = 1/9 and the maximum value at r = 1 
given again by Amax = In 3. 

For other values of the maximum Lyapunov exponent one 
can use a similar method starting from a piecewise linear map 
with the desired value of the Lyapunov exponent. On the other 
hand, if a constant Lyapunov exponent is needed, we can use 
instead the method leading to ( fT3] l but starting from an appro- 
priate chaotic map, such as ( |2TI ). 
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